State Tomography of a Qubit through Scattering of a Probe Qubit 
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We discuss the state tomography of a fixed qubit (a spin-f/2 target particle), which is in general 
in a mixed state, through f D scattering of a probe qubit off the target. Two strategies are presented, 
by making use of different degrees of freedom of the probe, spin and momentum. Remarkably, the 
spatial degree of freedom of the probe can be useful for the tomography of the qubit. 

PACS numbers: 03.65.Wj, 03.65.Nk, 72.10.-d 



I. INTRODUCTION 

The determination of a quantum state is a highly non- 
trivial problem. A wave function itself, or more generally 
a density operator, is not an observable and cannot be 
measured directly. One can see the state only through 
measurable quantities, which are related to the matrix 
elements of the density operator. It is therefore an inter- 
esting and important issue to discuss how to infer a given 
quantum state from a list of observed quantities. Such a 
problem is called state reconstruction or state tomogra- 
phy [1]. 

Recent advances in technology have been enabling us 
to engineer variety of peculiar quantum states to explore 
fundamental aspects of quantum mechanics. Generation 
of highly quantum states is one of the key elements for 
the realization of the ideas of quantum information [2]. 
It is clear that such issues cannot go without supports by 
the state reconstruction technology. Probing nontrivial 
correlations in a quantum system also helps us to clarify 
intrinsic characters of various many-body systems. 

From a practical point of view, direct accesses to tar- 
get systems are often limited. In such a case, we have 
to explore a way to probe the target in an indirect way. 
Scattering has always been considered a very powerful 
way to investigate many physical systems in a wide range 
of fields of physics, from elementary-particle physics to 
condensed-matter physics. Loosely speaking, all the 
physical processes to access targets can be regarded as 
scattering processes. 

In this paper, we focus on the state tomography of a 
qubit (spin-1/2 particle) via scattering, in a simple ID 
setup. We send a probe qubit to the fixed target qubit 
and see the state of the probe after its scattering off the 
target. Its scattering data contain the information on the 
target state before the scattering, from which we recon- 
struct the state. In order to reconstruct the qubit state, 
three independent scattering data are required. One pos- 
sibility is to make use of the spin degree of freedom of the 
probe. A collection of three transmission/reflection prob- 
abilities with three different sets of the initial and final 



spin states of the probe provides sufficient information 
for the reconstruction of the target qubit state. 

Notice here that in our setup, there is an additional 
(spatial) degree of freedom available in the scheme for 
the tomography of the spin, i.e., the momentum of the 
probe. It will be shown that this degree of freedom can 
be utilized to "optimize" (in the sense prescribed below) 
the scheme. It is interesting to see that the spatial degree 
of freedom can play a central role for the tomography of 
the spin. It will be demonstrated that three scattering 
data required for the qubit tomography are available by 
arranging different incident momenta and scattering di- 
rections, with the initial and final spin states of the probe 
being fixed. 



II. SETUP 

Suppose that a qubit A is fixed at x = on a ID line. 
We are going to discuss the reconstruction of its state, 
which is in general a mixed state pa, through the scatter- 
ing data of a probe qubit X off the target A. We assume 
that the following Hamiltonian describes the scattering 
process: 



p 

H = \-g(cr x ■ er A )S{x). 

2m 



(1) 



Here x and p{= —ihd/dx) are respectively the position 
and the momentum of the probe qubit X of mass m, and 
the potential produced by the fixed qubit A is assumed 
to be represented by the delta function. When X is scat- 
tered by A, their spins interact with each other through 
the Heisenberg-type interaction with a positive coupling 
constant g, with <jj being the Pauli operators acting on 
qubit J(= X, A). Sec Fig. 1. 

Let us start by solving the scattering problem of the 
Hamiltonian (1). The probe X is sent from the left (x < 
0) to the target A with a fixed incident wave number 
k (> 0), and scattered. The scattering (S) matrix clement 
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FIG. 1: A probe qubit X is sent with a fixed wave number 
k from the left to a target qubit A fixed at x — 0, scattered 
by the delta-shaped potential produced by the target qubit A 
with a spin-spin interaction of the Heisenberg type, and de- 
tected on the right or on the left with spin-sensitive detectors. 



is given by [3] 

(k'C\s\kO 

= 8(k' - k)5 cc - 2 m 8(E k , - E k )(k'('\V\* k 0, (2) 

where | kQ is the eigenstate of the free Hamiltonian Ho = 
p 2 /2m and \^kC) that of the total Hamiltonian H = H$ + 
V, both belonging to the same eigenvalue E k = h 2 k 2 /2m, 
with £ denoting the spin degrees of freedom of two qubits. 
The coordinate representation of the latter reads as [3] 



(a|**C) 



2m 



= (x\kQ - / dx' G k (x - x ')—V{x'){x'\y k Q, (3) 

where G k (x) = i(2fc) _1 e lfe ' x ' is the retarded Green func- 
tion in ID. It is easy to calculate the source term at 
x' = 0, 

<0|*feC) = <0|fcC) - itt{cr x ■ «r A )<0|tt fc O 

= (l + iCL<T X -cr A )- 1 {0\kQ, (4) 

with a dimcnsionlcss parameter 



n 



h 2 k' 



(5) 



Recall that the Heisenberg-type coupling cr x • ci is 
rewritten as 

ctx ■ ct a = -37>i + T 3 (6) 
in terms of the projection operators 

\ - a x ■ a a „ 3 + a x ■ <ta 

Vi i , 7>3 i , (7) 

onto the spin-singlet and -triplet eigenspaecs, respec- 
tively. This allows us to evaluate the inverse operator 
as 
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1 + iVLtJx ■ cta 1 — 3if2 1 + ifl 
We therefore obtain 

(k'C\S\K) = 8(k' - k)(C\T\C) + 8(k' + k)(C'\R\0, 



(8) 



where the scattering matrices responsible for transmis- 
sion and reflection, T and R, read 



R = nVi + r 3 V 3 
1 



(l-3iQ)(l + if2) 
T = 1 + R = t x Vx + t 3 V 3 
1 



(-3fJ 2 - itia x ■ <ta) 



(9) 



(i-3in ) (i + ^ 1 - 2in - <n ^-^ (10) 

Here the coefficients ru 3 ) and ti(3) coincide with the 
reflection and transmission amplitudes calculated sepa- 
rately for the spin-singlet(triplet) eigenspace, in which 
the interaction Hamiltonian is given by a scalar (i.e., not 
spin-dependent) potential —3g8(x) for the spin-singlet 
case or by g8(x) for the spin-triplet case. This implies 
that the problem can be reduced to an ordinary scatter- 
ing problem of a spinlcss particle. 



A. Scattering probabilities 

In the tomographic schemes we are going to discuss in 
the following, the probe X is sent with its spin polarized 
to the direction specified by a unit vector rii, and we see 
the probability of X being transmitted or reflected with 
its spin rotated to nj. Such probabilities are given by 

P = Tr{^/)x(V/|5(|^)x(^| ® Pa)&}, (11) 

where \ipi) x and \ipf) x are the incident and final spin 
states of X, which are expressed in the Bloch-sphere for- 
malism as 



\i>i)x(ipi\ = 



1 + rii ■ er x 



1 + rif ■ cr x 



(12) 

and S is to be substituted by T or R depending on 
whether X is transmitted or reflected. The target qubit 
A is in general in a mixed state pa, which is characterized 
by a Bloch vector v such that 

p A = 1 + V 2 <TA (0<H<1). (13) 

The transmission probability for a given set {rii, rif ,k} 
reads as 

p* = Of(n) + a ' t (n)(n f ■n i ) 

+ b t (Q)[3(rif ■ v) + (rii • v)] + c t (£l)v ■ (rif X rii) 

(14) 

with 

l + 7il 2 



at(n) = 

c4(Q) = 



2(l + fi 2 )(l + 9fi 2 )' 

i + zn 2 

2(1 + S1 2 )(1 

n 2 



(l + ft 2 )(l + 9ft 2 )' 
c t (Q) = -- 



(l + fi 2 )(l + 9ft 2 )' 



(15) 
(16) 
(17) 
(18) 
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These coefficients are plotted in Fig. 2 as functions of the 
incident wave number k of the probe qubit X. Similarly, 
the probability for the reflected case reads 

P r = a r (fl) + a' r (fi)(n r rii) 

+ b r (Q){(nf ■ v) — (n.i ■ v)] + c r (Q)v ■ (rif x rii) 

(19) 

with 



a r (fl) 
aj.(fi) 
b r {Q) 
c r (f2) 



3S1 2 (1 + 
2(l + fi 2 )(l + 9f7 2 )' 

» 2 (1 - 9Q 2 ) 
~2(l + n 2 ){l + 9fl 2 )' 

n 2 

(l + 2 )(l + 9^ 2 )' 
(l + 2 )(l + 9fi 2 )' 



(20) 
(21) 
(22) 
(23) 



which arc plotted in Fig. 3. Observe that the trans- 
mission and reflection probabilities P* and P r are both 
spherically symmetric in spin space. This reflects the 
symmetry of the Heisenberg coupling crx • &a- n o pre- 
ferred direction is present in the system. 



III. STRATEGIES FOR THE STATE 
TOMOGRAPHY 



Let us now discuss the tomography of the state pa, 
through the scattering data P t r . Our objective is to de- 
termine the three independent components of vector v, 
which exactly corresponds to the complete specification 
of the state pa, see (13). From an experimental point of 
view, this means that we need to arrange three indepen- 
dent experimental setups. As is clear from (14) and (19), 
the scattering probabilities P t r explicitly depend on the 
initial and detection orientations rii and rif of the spin of 
the probe qubit X as well as on the parameter Q contain- 
ing the incident wave number k and the coupling constant 
g; we can consider different strategics for the tomogra- 
phy, by properly choosing these parameters. In the next 
subsections, we discuss two of such tomographic strate- 
gies, in which the spatial and spin degrees of freedom of 
X play different roles. 

In the first approach presented in Sec. Ill A, a central 
role is played by the spin degree of freedom of X: we tune 
the incident and detection orientations rii and rif of its 
spin, while the incident wave number k is fixed. In order 
to completely reconstruct pa, we arrange the orientations 
rii and rif in three independent ways. For instance, we 
can choose three different detection orientations n f with 
the orientation of the incident spin rii being fixed, or we 
can tune both of them at the same time. The remaining 
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degree of freedom, the wave vector k, also plays an active 
role: we can "optimize" in the sense described below the 
tomographic scheme by tuning it appropriately. 

In Sec. IIIB, we discuss another approach, in which 
we fully make use of the spatial degree of freedom k. 
The tomography of the spin state pa is possible through 
the spatial degree of freedom. We collect three indepen- 
dent scattering data for different incident wave numbers k 
and scattering directions (i.e., transmission or reflection) 
when the incident and detection orientations of the spin 
rii and rif are fixed, from which wc reconstruct the spin 
state pa- Furthermore, one can follow the same criterion 
introduced for the first strategy in order to "optimize" 
the tomographic scheme. 



A. Strategy I: tuning m and n/ with k being fixed 

In this section, we reconstruct the initial state of the 
target qubit A by tuning the incident and detection ori- 
entations «.j and n/ of the spin of the probe qubit X 
with a fixed incident wave number k. In particular, we 
discuss two possible examples of this first tomographic 
strategy. 

Let us first fix the orientation of the spin of the inci- 
dent qubit X, m, and choose three different detection 
orientations rif = ni,n2,n3, which can be chosen to 
be orthogonal to each other. For instance, choosing ri\ 
perpendicular to rii, 



no 



rii x m, n 3 = rii 



(24) 



would be a natural choice as a reference frame. The 
transmission probabilities associated to these detection 
orientations with a fixed wave number k read 



3b t -c t b t \ hi 
ct 3b t h \ \v 2 
ibt) \v 3 









-1 








where 



(a = 1,2,3) 



(26) 



are the three independent components of the target vec- 
tor v along the axes of the reference system {rii, n 2 , 
introduced above and P l a the transmission probability 
when the spin of X is detected along direction n a . The 
three components of the vector v are readily obtained by 
inverting the relation (25) as 



(27) 









( a t 


M t 






a t 








\a t + a' t 



ith 



M t (0) 



l + O 2 
40 2 



'120 2 -40 0(1-30) 
40 120 2 -0(1 + 30) 
1 + 90 2 



(28) 



which completes the reconstruction of the state pa- Ob- 
serve that the only condition needed to be satisfied in 
order to invert the matrix in (25) to obtain Mt in (28) is 
< k < 00, which can be considered as a self-consistency 
condition for the present tomographic scheme. 

Another example of this first strategy, which is simpler 
from a computational point of view, is detecting the spin 
X oriented in the same direction as the incident spin, 
rif = rii . In this case, the transmission probability reads 



P 



v), 



+ 46 t (n, 

from which wc immediately obtain 

P'(l + 2 )(l + 90 2 ) - (1 + 50 2 ) 

Tlj • V = ttt; . 

40 2 



(29) 



(30) 



By choosing three different orientations for rii = w/, we 
gain the complete information on the vector v, that is, 
on the state pa- Observe that also in this case the self- 
consistency condition reads < k < 00. 

Similar schemes are available with reflection probabil- 
ities P r . Notice however that, if we choose rif — rii 
as we did in the above two schemes, no information on 
v is attainable from the reflection probability P r , see 
(19), where v disappears from P r for = rij. A pos- 
sible solution to this problem is to flip the orientations 
n a — > — n a (a = 1,2,3) in the choice of the reference 
system {ni,ri2 : n 3 ] in (24) for the former scheme, while 
measuring the reflected X in — rii direction instead of 
rii for the latter. In this way, the schemes with the re- 
flection probabilities work similarly to the ones with the 
transmission probabilities, with the same self-consistency 
condition needed for the tomographic reconstruction of 
the initial state of the target qubit, pa- 

The reason why these tomographic schemes that make 
use of the spin degree of freedom of X work is the fol- 
lowing. Let |t(|)) denote the state with spin parallel 
(anti-parallel) to the orientation rii of the incident probe 
spin X. Then, the state of A is in general expressed as 

PA = PnmA(1\+ Pllll) A(l\+P U mA(l\+ PU\1) A (% 

. (31) 

and the spin state of the probe X after the transmission 
by A becomes (apart from the normalization) 



Px 



«Tr A {r(|t)x(tl®PA) Tt} 



oc p TT + 



l + O 2 



jPll 



2iO 
+ 3zO 



h90 2 

puIOxUI 



ItMth 

2ifl 



40 2 



l-3iO 



1 + 90 2 

pitUMTI 



Pll\i)x(l\ 
(32) 



A similar expression is available for the reflection case. 
Observe that the component p^ of pa is associated to 
|t)x(tl m Px, Pn to lt)xU|, and so on: the spin state 
Pa is more or less "transferred" to X after the scattering. 
This is due to the Heisenberg coupling crx 1 &Ai which 
"swaps" the states between A and X. This is why we 
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can see the spin state of A by looking at the spin state 
of X. 

Until now, we have shown how the initial state of the 
target qubit A can be reconstructed by sending a probe 
qubit X with a fixed wave numebr k. A natural question 
would then arise: can we use the spatial degree of free- 
dom of X to optimize the above tomographic schemes in 
the sense that possible errors in the scattering data can 
least affect the determination of the state? It is actually 
possible. For instance, one would be able to reduce the 
effects of possible errors in the observations of the prob- 
abilities P t r on the reconstructed vector v, by properly 
tuning the incident wave number k. Observe that (27) 
and (30) arc linear mappings between the Bloch sphere 
of radius 1 and the probability space, associated to the 
scattering data P t,r . The reconstructed vector v is least 
sensitive to the errors in the observed probabilities P t r , 
when the volume of this probability space is maximum. 
Stated differently, under such a condition, the probabili- 
ties P* ,r are the most sensitive to the Bloch vector v to be 
reconstructed, and one can do a better tomography. For 
the first scheme, in which we tune the detection orienta- 
tion rif of the spin of the probe X, the volume associated 
to the probabilities is maximum when the Jacobian of the 
map (27), which is given by the determinant 



2\2 



detAf t (fi) = 



(i + n a ) d (i + 9n a ) 

40^ 



(33) 



is minimum. For the second scheme, in which the inci- 
dent and detection orientations of the probe spin X are 
the same rif = rij, the coefficient multiplying the scat- 
tering probability P* in (30), 



A(n) 



(l + » 2 )(l + 9f2 2 ) 
4^ 2 



(34) 



is to be minimized. In Fig. 4, [det M^fi)] 1 / 3 and X(fl) are 
plotted as functions of the incident wave number k of the 
probe X. They become minimum at 1/ft = h 2 k/mg — 
yl + \/217 /2 ~ 1.98 and y/3, respectively. In particu- 
lar, observe that, at h 2 k/mg = -\/3, the formula for the 
tomography (30) is reduced to 



n, ■ v 



4P* - 2, 



(35) 



from which it immediately follows that the probability 
p t 



ranges 



I<P'<3. 
4 - - 4' 



(36) 



which is the maximum in this scheme. The same anal- 
ysis can be applied to the schemes with reflection prob- 
abilities, for which the optimal momenta are h 2 k/mg = 

J (—3 + \/33)/2 — 1.17 and \/3 for the two schemes, re- 
spectively. 
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[det Mtin)] 1 /' 3 



h 2 k/mg 



FIG. 4: The Jacobian of the linear map (27) between the 
probability space and the Bloch sphere, det Mt(fi) given in 
(33) (solid), and the coefficient A(f2) (dashed), as functions 
of the incident wave number k of the probe X, which are 
both to be minimized by tuning k for the optimizations 
of the two approaches of Strategy I, presented in the text. 
[detM t (fi)] 1/3 is actually plotted instead of detM t (fi), to 
better compare it with A(f2). [det Mt(f2)] 1//3 is minimum at 
h 2 k/mg — y] 
dots). 



T27/2, while A(fi) at y/3 (indicated by 



B. Strategy II: tuning k with n< and rif being fixed 

The spatial degree of freedom of X can itself play a 
fundamental role for the tomographic reconstruction of 
the target spin state pa- We fix the incident and detec- 
tion orientations rii and rif of the qubit X and tune the 
wave number k to collect sufficient number of scattering 
data required for the tomography of pa- 

Let us fix rif perpendicular to rii and select three dif- 
ferent incident wave numbers k = ki, k 2 , k 3 . A collection 
of the three transmission probabilities with these wave 
numbers yields 



'P'(tti)' 
PW 



'3&t(ni) bt(Sli) ct(fii)\ fv f 

3&t(n 2 ) b t (fl 2 ) ct(fia) Vi 

v 36 t (0 3 ) b t (fl 3 ) c t (n 3 )J \v± 

'ot(ni)\ 

a t (« 2 ) , (37) 

Mas)/ 



where Q,j = mg/h 2 kj (j = 1, 2, 3), and 

Vf = rif ■ v, Vi = rii - v, v± = (rif x m) ■ v (38) 

are the three components of the vector v with respect 
to the reference frame fixed by {nf,rii,nf x n^}. This 
relation cannot be inverted, since two columns of the ma- 
trix are proportional to each other, irrespectively of the 
choice of {^1,^2,^3}. This is due to the fact that, once 
the orientations rii and rif are fixed, only two compo- 
nents of the vector v along the directions specified by 
3n f + rii and rif x rii are involved in the transmission 
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probability (14), while it is insensitive to the other com- 
ponent of v perpendicular to the plane spanned by these 
directions. The same happens if we collect the reflection 
probabilities (19) in a similar way. This problem how- 
ever can be overcome by combing the transmission and 
reflection probabilities. For instance, 



,i*(fi2), 



'& r (ftl) -& r (fii) 

36 t (fii) 6t(fii) 
v 3& t (0 2 ) bt(Sl 2 ) 

a r (r^i > 



Cr(^l) 
C t (fll 
Ct(fl 2 




(39) 



which is inverted as 

V 

»i | =JV(ni,n a ) 



^p r («i) 
p'(^i) 



V(fii) N 

a t (fii) 
,a t (ft 2 F 



(40) 



with an inverse matrix iV(f2i,fi 2 ), when 



(i + r> 2 ) 2 (i + 9r> 2 ) 2 (i + np(i + gn|) 
4fif^ 2 (fii -n 2 ) 



detiV(fti,fi 2 ) = - 

(41) 

is nonvanishing and finite. This condition is fulfilled by 
fci and fc 2 that are both finite and different from zero as 
well as different from each other. In this way, we can re- 
construct the spin state pa by making use of the spatial 
degree of freedom of X. Note that only two experimental 
setups are actually needed to collect the three scattering 
data: the two probabilities P*(f2i) and P r (rii) are ob- 
tained at the same time in a single setup, by sending the 




ft ki/mg 

FIG. 5: (Color online) Contour plot of |det JV(Oi, 2 )| 1/3 
given in (41), as a function of k\ and /c2. Only the 
contours in the range |det N(Qi, O2) | 1 / 3 < 20 are shown. 
|det N(fli, f^)! 1 / 3 takes its minimum value ~ 8.74 at 
(Jfei, fa) d (1.51, 5.13) in the unit mg/fi 2 (indicated by a dot). 



probe X with ki and seeing whether it is transmitted or 
reflected with its spin orientated to nj. 

Similarly to the first strategy discussed in the previous 
subsection, the present scheme is optimized by appropri- 
ately choosing the two incident wave numbers k\ and fc 2 
of the probe X. Under the same criterion as the one for 
the first strategy [maximizing the volume of the proba- 
bility space associated by the linear mapping (39) to the 
Bloch sphere of radius 1 , to which the vector v belongs] , 
the optimal choice of (fci, fc 2 ) is found by minimizing the 
quantity |det A^(Sli, S^ 2 )| 1/3 from (41). See Fig. 5, where 
the optimal set of wave numbers for the present strategy 
is found at (fcx, fc 2 ) — (1.51, 5.13) in the unit mg/h 2 . 



IV. CONCLUSIONS 



In this paper, we have discussed the state reconstruc- 
tion/tomography of a fixed qubit through the scattering 
data of a probe qubit off the target. We have presented 
two different strategies for the tomography, in which the 
spin and spatial degrees of freedom of the probe qubit 
play different roles. 

The first strategy makes use of the spin degree of free- 
dom of the probe. The spin state of the target is more or 
less transferred to the probe spin during the scattering, 
and therefore, we can infer the spin state of the target 
through the state tomography of the probe spin. The 
other degree of freedom, the momentum of the probe, 
can be utilized to optimize the tomographic scheme. 

The spatial degree of freedom can also play a central 
role for the state tomography of the target spin. In the 
second strategy, three scattering data required for the 
state tomography of a target spin are collected by choos- 
ing different incident wave numbers and scattering di- 
rections (transmitted or reflected), with the incident and 
detection orientations of the probe spin being fixed. This 
tomographic scheme can be optimized also by appropri- 
ately tuning the set of the incident wave numbers. 

The strategies introduced in this paper for a single 
fixed qubit can be generalized to the tomography of mul- 
tiple qubits. In particular, the detection of entanglement 
would be an important task in the light of quantum infor- 
mation. In order to reconstruct the state of N spins, we 
need 4 N — 1 different experimental setups. Imagine, for 
instance, how to choose 15 different sets of orientations of 
the incident and detection spin states of the probe qubit 
for the two-qubit tomography. If the number of the tar- 
get qubits grows, tuning the probe spin to the different 
orientations required for the first strategy would become 
more and more difficult. In such a case, the second strat- 
egy may provide a way out of this problem. It would be 
worth exploring such a potential of the scheme, which 
would be an interesting future subject. 
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